Model analysis of thermal UV-cutoff effects on the chiral critical surface 
at finite temperature and chemical potential 
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We study the effects of temporal UV-cutoff on the chiral critical surface in hot and dense QCD 
using a chiral effective model. Recent lattice QCD simulations indicate that the curvature of the 
critical surface might change toward the direction in which the first order phase transition becomes 
stronger on increasing the number of lattice sites. To investigate this effect on the critical surface in 
an effective model approach, we use the Nambu-Jona-Lasinio model with finite Matsubara frequency 
summation. We find that qualitative feature of the critical surface does not alter appreciably as we 
decrease the summation number, which is unlike the case what is observed in the recent lattice QCD 
studies. This may either suggest the dependence of chemical potential on the coupling strength or 
due to some additional interacting terms such as vector interactions which could play an important 
role at finite density. 
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I. INTRODUCTION 

The study of the chiral critical point (CP) in the 
phase diagram of hot and dense quark matter is one of the 
central issues in Quantum Choromodynamics (QCD) [l[ . 
While, it is widely accepted that the QCD phase tran- 
sitions concerning chiral symmetry restoration and color 
deconfinement are crossovers with increasing tempera- 
ture T for small chemical potential fi ~ 0, the order of 
the phase transitions along the /j, direction for small T 
is still under considerable speculation. Model analysis 
indicate that a first-order phase transition occurs with 
increasing fj, and for small T [2[ . The above observations 
lead us to expect the existence of a critical point located 
at the end of the first order line in the phase diagram for 
some intermediate values, Te and he- 

A first principle determination of the phase diagram 
by solving QCD itself is difficult due to the strongly in- 
teracting nature of matter at low-energies/temperatures 
which significantly restricts the range of applicability of 
perturbative calculations. We must then rely on non- 
perturbative techniques such as lattice QCD (LQCD) or 
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some low-energy effective theories of QCD. The LQCD 
simulations are known to be a viable approach for micro- 
scopic calculations in QCD, and have recently reached a 
reliable level at finite T and fi = [3]. However, these 
simulations are not yet able to provide a conclusive un- 
derstanding of the QCD phase diagram due to severe 
limitations posed by the well-known "sign-problem" at 
finite /x, and difficulties dealing with small quark masses, 
though only very approximate methods are available for 
simulations at small /i values fi]. Thus, it is important 
to develop low-energy effective models that may show 
consistency with lattice results and can be extrapolated 
into regions not accessible through simulations. Among 
them, the local Nambu-Jona-Lasinio (NJL) model Q 
and its proposed extended version to include coupling 
of quarks to Polyakov-loops, the so-called Polyakov-loop 
NJL (PNJL) model @ , are useful to study the quark sys- 
tem at finite T and/or fi. Such effective models share the 
same symmetry properties of QCD and successfully de- 
scribe the observed meson properties and chiral dynamics 
at low-energies (see, e.g., |7Hlfj|V 

Using the framework of (P)NJL model we search the 
CP and analyze the order of the chiral phase transition 
by varying the current quark masses (m u , md, m s ). We 
usually set, for simplicity, to<j = m u and investigate the 
phase transitions in the m u -m s -fj, space. Renormaliza- 
tion Group (RG) analysis of chiral models at /i ~ con- 
clude that there is no stable infra-red (IR) fixed-point 
for quark flavors Np > \/3 [ff| , indicating that the ther- 
mal phase transition is of fluctuation induced first or- 
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der for two or more flavors realized in the chiral limit 
w«,rf.s = 0. However, the transition becomes a crossover 
for intermediate quark masses because of explicitly bro- 
ken chiral and center symmetries. Hence, it is naturally 
expected that there should be a "critical boundary" sep- 
arating the regions of the first order phase transition and 
crossover between small and intermediate quark masses. 
Although LQCD results support the above model pic- 
ture qualitatively, there had been a huge quantitative 
difference between the two kinds of analyses, even at 
zero chemical potential where the LQCD does not suf- 
fer from the sign-problem; the value of the critical mass 
obtained in the (P)NJL model is about one order of mag- 
nitude smaller than the value in the LQCD analyses. In- 
spired by recent works reporting that the critical mass 
may become smaller when the number of lattice sites is 
increased [IH, [l]| , the present authors studied the criti- 
cal boundary in the (P)NJL model with finite Matsubara 
frequency summation N at zero chemical potential [l4j |. 
There it was found that the critical mass actually be- 
comes larger if one decreases the Matsubara summation 
number N (see, Eq.Q), thereby showing the correct ten- 
dency to explain the quantitative difference between the 
LQCD and (P)NJL model results. 

In this Letter, our goal is to study the critical bound- 
ary for all values of the chemical potential, which may 
eventually tell us the location of the CP in the QCD 
phase diagram. As already mentioned that since an ab 
initio determination of the critical point in QCD is a dis- 
tant hope, nevertheless, it is worth studying the (P)NJL 
model with finite frequency summation at finite chem- 
ical potential that should capture the essential qualita- 
tive features of the results expected in lattice studies. 
Note that in the m u -m s -[i space, the critical boundary 
becomes a surface called the "critical surface" . More pre- 
cisely, through this study we would like to investigate the 
qualitative behavior of this critical surface whose shape 
can critically determine whether the CP exists. We can 
then compare our results with the recent lattice predic- 
tions. 



fermion interaction term called the Kobayashi-Maskawa- 
t'Hooft interaction whose coupling strength is <7d 0- 
The subscripts (i, j) indicate the flavor indices and the 
determinant runs over the flavor space. This term is in- 
troduced to explicitly break the Ua(1) symmetry. 

To study the chiral dynamics, we solve the gap equa- 
tions which are derived through differentiating the ther- 
modynamic potential 51 by the order parameters of the 
model: 



on 

dm* 







u, d, s 



(4) 



The order parameters m* are the constituent quark 
masses. Since, we set nid — m u in our analysis, this 
should lead to the isospin symmetric result m* d = m*, 
reducing the number of gap equations from three to 
two. The thermodynamic potential Vt is defined by 
£1 = —\nZ/(f3V), where Z is the partition function, 
/3(= 1/T) is the inverse temperature and V is the vol- 
ume of the thermal system. 

In the mean-field approximation, after some algebra, 
we arrive at the following expressions for the gap equa- 
tion: 



m* u = m u + 2ig s N c trS u - 2g B JV c 2 (trS u )(trS s ), 
m* =m s + 2ig s N c trS s - 2 ffD A 2 (trS m ) 2 , 



(5) 



where N c (= 3) is the number of colors and tvS 1 is the 
chiral condensate written explicitly as 



i ttS* = 4m* 



d3p (IT) ^ 



(2tt) 
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Here w n = 



7rT(2n + 1) is the Matsubara frequency and 
+ m* 2 is the energy of the quasi-particle. A 



detailed calculation for deriving the gap equations Eq.© 
is clearly presented in the review paper [8[ . 

To study the UV-cutoff effects in the model, we cut 
the higher frequency modes in the Matsubara sum as 
employed in |14| . 
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II. MODEL SET UP 

The NJL model Lagrangian in the 3 flavor system is 
written by 



£njl 



q (i$ - rh) q + Ci + C 6 , 

8 

9s 
2 



a=0 



£e = -go [dct^(l - j 5 )qj + h.c. ] . 



(1) 
(2) 
(3) 



Here rh is the diagonal mass matrix (m„, m^, m s ) in the 
flavor space which explicitly breaks the chiral symmetry. 
£4 is the 4-fermion contact interacting term with cou- 
pling constant , and A a is the Gell-Mann matrix in the 
flavor space with Aq = y / 2/ 7 3diag(l, 1, 1). £g is a 6- 



E - E 



(7) 
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This model has five free parameters {m u , m s , A, gs, gr>}' 
two current quark masses, a 3-dimensional momentum 
cutoff, a four-fermion and a six-fermion coupling con- 
stants. Following Q, we set m u — 5.5MeV fixed for all 
values of N, while the remaining four parameters are fit- 
ted each time by using the following physical observables 



m OT = 138 MeV, 
m K = 495.7 MeV, 



U = 93 MeV, 



nty = 958 MeV. 



The parameter fitting for various N has been done in [lj] , 
and we employ the same values in this analysis which 
are again displayed in Tab. U for the convenience of the 
reader. 
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N 


m 3 (MeV) 


A (MeV) 




ff D A b 


15 


134.7 


631.4 


4.16 


12.51 


20 


135.0 


631.4 


4.02 


11.56 


50 


135.3 


631.4 


3.82 


10.14 


100 


135.4 


631.4 


3.75 


9.69 


coo 


135.7 


631.4 


3.67 


9.29 



TABLE I: The various fitted parameters for different N [jj 



III. CHIRAL CRITICAL SURFACE 

The main purpose of this Letter is to determine 
the chiral critical surface in the (P)NJL model with fi- 
nite Matsubara summation. The critical surface is the 
set of all critical points in the m u -m s -fi space which are 
analyzed by scanning the space for discontinuities of the 
chiral condensate. It should be noted that for each value 
of N, we treat both the current quark masses m u and 
m s as free parameters in obtaining the critical surface 
once the other parameters, namely, A, <?g and <?d are de- 
termined by fitting to the physical parameters, as shown 
in Tab. Q] Of course, we will eventually be interested 
in the case of the real (physical) current quark masses 
m u ~ 5.5MeV and m s ~ 136MeV, in order to determine 
the possible existence/non-existence of the CP through 
our model analysis. Because our motivation is to make a 
direct comparison of our results with that of LQCD where 
the simulations are mainly performed in the m u = m s 
symmetric case at finite /i, we shall also consider this 
case. In the actual numerical calculations, we scouted 
out the critical masses (m uc ,m sc ) for each fx by searching 
for discontinuities in the solutions of the gap equations 
Eq.([5]) in the entire m u -m s -fj, space. 

The LQCD and model studies indicate a crossover re- 
alized at /i = for physical current quark masses. This 
means that the curvature of the critical surface will tell 
us whether the CP is favored in the phase diagram. To 
be more concrete, if the region of the first order phase 
transition expands with increasing /x, the physical quark 
mass line will intersect with the critical surface and this 
will end up as a CP. If on the other hand, the first or- 
der phase transition region shrinks with /i, there is less 
chance of an appearance of a the CP and a crossover 
transition will be favoured for the whole range of T and 

fl. 

In the LQCD calculations, the curvature of the critical 
surface along the m u = m s symmetric line is analyzed 
by obtaining the critical mass m c through the following 
Taylor expansion formula 



m c (0) 



1 



fe=i 



2 A- 



(8) 



which so far yielded the following results: c\ = —3.3(3), 
c 2 = -47(20) for N t = 4, and c x = 7(14), -17(18) (pre- 
liminary) for a leading order (LO) and next-to-leading or- 
der (NLO) extrapolation in fi 2 , respectively, for N t — 6 
where Nt represents the number of the lattice sites in 



the temporal direction [l3(. These results are graphi- 
cally represented in Fig. [3] and [5] in the following sections 
which indicate that the sign of the curvature has not yet 
been determined from lattice simulations. 
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IV. NJL MODEL RESULTS 



The phase diagrams of the NJL model for N = 15, 
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FIG. 1: Position of the CP in the T-^t phase diagrams from 
the NJL model with finite Matsubara summation for N = 
15, 50, oo. The dotted and the solid lines represent crossovers 
and first order phase transitions, respectively. The N = oo 
case corresponds to the traditional NJL model. 

50 and oo, respectively, are shown in Fig. [TJ In fact, 
we have also studied the model for other values of N, 
however, there were no significant qualitative differences 
and we prefer to select just the above three representa- 
tive cases for graphical clarity. Here, it may be notewor- 
thy mentioning that for values of TV < 15, it becomes 
a matter of numerical challenge to perform simulations 
to determine the phase boundaries. So, henceforth, we 
shall be displaying our results only for the above three 
cases. Note that the case N = oo corresponds to the 
traditional NJL model. In drawing the phase diagrams, 
we apply the same criterion employed in [To| where the 
phase transition or crossover is defined by the condition, 



(uu) 



/To 



1 

2' 



(9) 



(uu)t being the expectation value of the chiral conden- 
sate for the up quark at temperature To and /i = 0. We 
choose To = for the N = oo traditional model, while 
we set To = 50MeV for the case with finite N since the 
model is ill-defined for small T, as discussed in [T3 |. This 
is why we choose not to display the results for the small T 
region where the curves are no longer physically reliable. 
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Here it is seen that the region below each of the curves 
that represents the chiral symmetry broken phase ex- 
pands with decreasing N. This comes from the fact that 
the coupling constants become larger with decreasing N, 
being consistent with RG arguments that the coupling 
strength becomes smaller when one considers the physics 
at higher momenta, i.e., for larger N. When the cou- 
pling strength grows the chiral condensate tends to en- 
large, which can be easily seen from Eq.©. Thus, it 
is naturally understood why the transition temperature 
and chemical potential increase with a smaller choice of 
the summation number N. 
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FIG. 2: m e (/Lt)/m c (0) vs /j,/(ttT) plot for different N values 
from the NJL model, along with the corresponding results 
obtained from lattice simulations for Nt = 4 and Nt = 6. 
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In Fig. [2j we show the numerical results for 
m c (/x)/m c (0) as a function of ^l/{ttT), along with the 
corresponding results obtained in the recent LQCD sim- 
ulations for Nt = 4 and 6 [IH, respectively. We see that 
the slope of the curves tends to go up only very slightly 
on decreasing N from N = oo down to iV = 15. On the 
other hand, the results from lattice simulations do not 
yield conclusive results so far, as clearly revealed from 
the above figure, where the sign of the curvature has not 
yet been constrained for the N t = 6 case. Thus, in com- 
parison with the lattice studies, we find that the quali- 
tative behavior of our model results hardly changes by 
using different N values in our analysis. This means that 
there is indeed a stark quantitative contrast between the 
NJL model results and the lattice predictions. 

As the final remark in this section, we display our re- 
sults for the critical surface obtained for the N = oo and 
15 cases in Fig.[3j It is clearly seen that, although the val- 
ues of the critical mass changes about factor of two, the 
qualitative picture docs not differ as wc change the value 
of N. Thus, the region of the first order phase transi- 
tion expands with respect to \x for all TV values, since the 
effect of variation of N is rather too nominal to change 
the sign of the curvature of the critical surface. This 



FIG. 3: The chiral critical phase surfaces for (a) N = oo, and 
(b) N = 15 from the NJL model in the m u -m s -fi space. 



would then mean that the NJL model will always favour 
the CP with physical current quark masses at some finite 
value of fj,, e.g., we obtain the CP when fi goes up around 
fi c = 324 (342) MeV for N = oo (15), as exhibited in the 
phase diagram Fig. [1] 



V. THE PNJL MODEL EXTENSION 

It is also intriguing to study the critical surface in the 
PNJL model, because of the closer resemblance to QCD 
which treats the chiral and deconfinement phase tran- 
sitions simultaneously. In the PNJL model, the order 
parameter for the deconfinement phase transition is the 
Polyakov-loop and it is described by a global mean-field 
being similar to the chiral condensation in the traditional 
NJL model. The Lagrangian of the PNJL model is writ- 
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ten by 

£pnjl = A) + £ 4 + U + $*, T) , 
C = q (i$ - i74 At -m)q, 
W($,$*,T) = -&T{54e" a/T $$* 
+ ln[l-6$$*+4($ 3 + $* 3 



(10) 



(11) 

3($$*) 2 ]}, (12) 



where W is the Polyakov-loop effective potential and $ 
and $* are the traced Polyakov- and the anti-Polyakov- 
loop, respectively. They are defined by $ = (l/N c )trL, 
$* = (l/JV c )trit with L = 7>exp[i ^AtAA and A 4 = 
LA . There are several candidates for the Poly akov-loop 
potential in defining the PNJL model |6j, |l(J, [l6|, and 
we adopt the strong-coupling inspired form of Eq.([T2|) 
following [!(}. In the above expressions, the parameter 
a solely parametrizes the strength of the Polyakov-loop 
condensate for the deconfinement phase transition, while 
the parameter b controls the relative strength of the mix- 
ing between the Polyakov-loop and chiral condensates, 
with a smaller value of b signifying chiral phase transition 
dominating over deconfinement. Here, the parameters a 
and b are set as a = 664 MeV, and &-A -3 = 0.03. Re- 
garding the rest of the model parameters, it is legitimate 
to use the same ones fixed in the NJL model because the 
Polyakov-loop extension is likely to affect the system only 
at finite temperatures comparable to the critical temper- 
ature T c . At much lower temperatures the chiral con- 
densate is only very marginally modified by the Polyakov 
loops. 



VI. PNJL MODEL RESULTS 

Let us now discuss the results in the PNJL model 
with finite frequency summation. 

In Fig. 21 we present the phase diagrams resulting from 
the PNJL model with N = 15, 50 and 00. Here we see 
that the curves are very similar to the ones in Fig. [TJ 
however, the critical temperatures are about a factor 
of two or more larger than those obtained via the NJL 
model. This is almost the same quantitative difference 
what is observed between the traditional (N = 00) NJL 
and PNJL models. 

We also display the corresponding curves for the ratio 
m c (/x)/m c (0) and compare the results with the LQCD 
simulations. Again, the nature of the results exhibit sim- 
ilar qualitative characteristics with the ones obtained in 
the NJL model; the ratio does not change appreciably. 
However, the ratios are numerically larger than that ob- 
tained with the NJL case. This result can be interpreted 
as an effect of the Polyakov-loops tending to suppress 
unphysical quark excitations below T c [Io| . 

Finally, in Fig. |H1 the critical surfaces in the PNJL 
model with N — 00 and 15 are displayed. We confirm 
almost the same qualitative features in these two figures, 
quite similar to that previously found in the NJL case. 
It is interesting to note that such qualitative similarity 
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FIG. 4: Position of the CP in the T-[i phase diagrams from 
the PNJL model with finite Matsubara summation for iV = 
15, 50, 00. The dotted and the solid lines represent crossovers 
and first order phase transitions, respectively. The N = 00 
case corresponds to the traditional PNJL model. 
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FIG. 5: m c (fi)/m c (0) vs fi/(irT) plot for different N values 
from the PNJL model, along with the corresponding results 
obtained from lattice simulations for Nt = 4 and Nt = 6. 



between the NJL and PNJL model results is a priori non- 
trivial, since the deconfinment phase transition order pa- 
rameters, namely the Polyakov-loops <I> and respec- 
tively, may cause the system with of chiral and decon- 
finment phase transitions to deviate significantly (both 
qualitatively and quantitatively) from a system with only 
chiral phase transition, especially in the vicinity of the 
CP. 
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FIG. 6: The chiral critical phase surfaces for (a) JV = oo, and 
(b) N = 15 from the PNJL model in the m u -m 3 -fi space. 

VII. SUMMARY AND DISCUSSION 

In this Letter, we studied the UV-cutoff effects on 
the chiral critical surface with two light and one heavy 
flavors using the (P)NJL model and found that its curva- 
ture is not appreciably affected by the finite Matsubara 
frequency summation number. On the other hand, the 



current lattice simulations are not decisive at the mo- 
ment which are beset with larger lattice cut-off effects 
than finite density effects, making continuum extrapola- 
tions doubtful. However, there is room for further pre- 
cise lattice calculations in future with greater number of 
lattice sites and development of new techniques for finite- 
density simulations that may be necessary to make more 
definite conclusions. 

As a final note, we point out the possibility of the 
temperature and density dependence of the coupling con- 
stants that may become important at high-energies, as we 
expect the couplings to run with respect to the energy 
scale. However, in this Letter, we used constant values 
of the couplings gs and gr> which were fixed from fitting 
to physical quantities once and for all at small temper- 
ature and zero chemical potential, to test the effects of 
the temporal UV-cutoff. These dependencies may turn 
out to have crucial effect on the critical surface when one 
considers the system at high densities, where it is ex- 
pected to be dominated by non-hadronic states. Thus, 
the location of the CP is indeed sensitive to the nature 
and magnitude of the coupling constants. In fact, it was 
actually found in that the Ua(X) anomaly strength 
modeled through the chemical potential dependent cou- 
pling go may change the curvature of the critical surface, 
as well as its sign, resulting in a characteristic "back- 
bending" of the critical surface as a function of pk. This 
reflects the fact that the density dependence of the cou- 
pling strengths plays crucial role when investigating the 
chiral critical surface. 
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